Depinning of elastic manifolds 
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We compute roughness exponents of elastic d-dimensional manifolds in (d + l)-dimensional em- 
bedding spaces at the depinning transition for d = 1, . . . ,4. Our numerical method is rigorously 
based on a Hamiltonian formulation; it allows to determine the critical manifold in finite samples 
for an arbitrary convex elastic energy. For a harmonic elastic energy (A^-model), we find values 
of the roughness exponent between the one-loop and the two-loop functional renormalization group 
result, in good agreement with earlier cellular automata simulations. We find that the A^-model is 
unstable with respect both to slight stiffening and to weakening of the elastic potential. Anharmonic 
corrections to the elastic energy allow us to obtain the critical exponents of the quenched KPZ class. 



Elastic manifolds in random media are an important 
issue of current research in statistical physics. H'S In the 
zero-temperature motion of these manifolds, subject to 
a driving force /, the "depinning threshold" /c plays a 
central role: For forces f > fc, the elastic manifold moves 
with finite velocity, while it is pinned for f < fc- Among 
the subjects studied at and around fc are the "creep" 
motion for / < at finite temperature, the scaling of 
the velocity for small forces f ^ fc, and the statistical 
properties of the pinned critical manifold at fc, especially 
its roughness exponent 

If one neglects velocity-dependent terms in the equa- 
tions of motion of the manifold, which one assumes to be 
a single- valued function /i(x,i) (no overhangs), its dy- 
namics is governed by a functional x}) incorpo- 
rating potential energy due to the driving force /, the 
disorder /^(x, h), as well as its internal elastic energy E^i, 

Note that x is a d-dimensional vector, in an embedding 
space of dimension d + 1. Dimensional analysis suggests 
that the harmonic approximation for Ed provides the 
only relevant term for the interplay between disorder and 
elasticity. This yields the quenched Edwards- Wilkinson 
(EW) equation 

dth{x,t)^f + r]ix,h) + aV^h{x). (2) 

The model described by eq. (||) has been an important 
testing ground for the concepts and techniques originally 
developed in the field of critical pJaenomena. Functional 
renormalization group tcchniquegaa were used to pertur- 
batively compute critical exponents of the quenched EW 
equation, which were believed to describe generic driven 
manifolds at the depinning threshold. 

According to the framework provided by the functional 
renormalization group, the manifold is flat = 0) for 
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FIG. 1: Schematic description of a d = 2-dimensional mani- 
fold, driven by a force density /. The (d -I- l)-dimensional 
embedding space is periodic in the manifold's dimensions 
ii, . ■ . ,id, but also in the variable h. 

d > due — 4. Below this upper critical dimension, the 
roughness exponent is expressed in aii,fe-expansion in e = 
4 — d. The first order (one- loop) termHu of the expansion, 
gives C = e/3. Initially, this result was Jpelieved exactcl 
for all d = 1, 2, 3. However, Chauve et alB have obtained 
the two-loop corrections, and found them to be non-zero. 

In this work, we compute the roughness exponent of 
critical manifolds in finite (d -I- l)-dimensional samples 
for d = 1,... ,4. Our approach is rigorously based on 
a Hamiltonian formulation, as we study directly the dy- 
namics of eq. (|l|) for a general convex elastic energy E^: 
A powerful numerical algorithm allows to solve this prob- 
lem for a slightly modified version of eq. (|l|), with a 
discretized vector x (x ^ i) and a continuous variable 
h. The short-range elastic energy depends on the next- 
neighbor distances Ai^^ between lattice point i and its 2d 
neighbors izt Sj, as shown in Fig. ^ 

^EE^'=i(i'^i-'^i+^.i)- (3) 

J = l ±Sj 

We introduce periodic boundary conditions for x and 
h (c/ Fig. mE Continuous, periodic random potentials 
are const*U|Cted from M normally distributed random 
variables .LfH In such a system, the sample-dependent crit- 
ical force fc is well defined: pinned configuration with 
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Wi = Vi do (do not) exist for driving forces smaller 
(larger) than f^. Our algorithnJa allows to decide quickly 
whether pinned configurations exist at a given driving 
force, and then zooms in on the critical force and the 
critical manifold. 

For the first time, we are thus able to unambiguously 
compute the object of most direct theoretical interest in 
this problem, namely th£ critical manifold, in dimensions 
d = 2, 3, 4. We computaa the roughness exponent ( from 
the disorder-averaged mean square deviations of the crit- 
ical manifolds h'^: 
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(4) 



As the width of the manifold is of order ^ L'^, we scale 
the lateral extension of the sample at least as M ~ L''. 

Especially in higher dimensions, direct simulations of 
eq. (^ are not viable, because they violate a crucial no- 
passing theorem JlJ which we respect. 

Because of the expected universality and the difficulties 
of direct simulation, attempts to compute the critical ex- 
ponents of the quenched EW equation considered discrete 
systems with non-har nili . oj iian dynamic rules. These cel- 
lular automata modclsMtJ arc not clearly conmpcted with 
an equation of motion in the continuum limittJ and have 
often given only rough estimates for ^. 

Nevertheless, simulations using cellular automata illus- 
trated that critical elastic manifolds in random media are 
not generally described by the quenched EW equatkm, 
but may fall into two broad classes: The automataB'til 
summarized in Ref. |l^ yield critical exponents which 
are close to those obtained by anaJbitical work on eq. (||); 
on the other hand the automataliJ^t3 summarized in Ref. 
^give very different values for the exponents, incompat- 
ible with the results obtained by. ren rpaialization group 
results, but closer to experiment .Il30't§ There has been 
much confusion about the actual values of these expo- 
nents for d > 1. It is not clear whether the upper critical 
dimension of this class is also d^c — 4. 

In our controlled Hamiltonian approach, we recover 
this very rich behavior as a dependence of the roughness 
exponent ^ on the functional form of the elastic energy: 
We investigate elastic energies of the form E^i oc |A|^, 
and refer to them as A'^'-models. The A^-model stud- 
ies directly the eq. (^. Models with 7 7^ 2 differ in 
their critical behavior from the harmonic model. A spe- 
cial role is played by the A'^-model, which represents the 
first non-harmonic corrections of a general clastic energy 
Sei ~ aAV2 6A4/12 . . . 

We first discuss this important issue for an elastic 
string {d = 1) in a 2-dimensional medium, which has 
been abundantly studied in the past. As shown in Fig. 
we find for the A^-model a value of the roughness ex- 
ponent = 1.26 ± 0.01, which compares well with the 
two- loop calculation of Chauve et ai, while clearly ex- 
cluding the one-loop result Ca2°'°°'' = 1- ^ value in ex- 
cess of one for the (1 -I- l)-dimensional roughness has 
been found in direct numerical integrations of the equa- 




FIG. 2: Mean square extension vs system size for the 
(1 + l)-dimensional elastic string. Upper curve: harmonic 
elastic energy _Eei(A) — M ~ L^'^. Lower curve: elas- 

tic energy Eei{A) — M = L . The inset shows the esti- 

mated value of C, for all the data with L > Lmin as a function 
of I/min (A^-model). The absence of systematic trends leads 
us to conclude that Ca^ = 1-26 ± 0.01, Ca-i = 0.635 ± 0.005. 



tions of motionjll in Monte Carlo simulations J§@ and in 
the cellular automata of Ref. |l^. 

Figure || also shows our results for the A''-model. We 
findEJ = 0.635 ±0.005. This exponent coincides with 
the value found in Ref. |l^. In lattice models, we already 
showed that the roughness exponent does not change for 
the A^-model and even in the presence of a metric con- 
straint (bounded |A|). In this sense, the A*-model is 
only one representative of systems with a stronger than 
harmonic elastic energy. The values obtained from Fig. |^ 
are much more precise than the earlier ones. 

In this problem, the study of rotational invariance of 
the equations of motion has occupied a crucial role. In 
fact, eq. (^, and therefore the equation of motion of our 
A^-model, is invariant under a tilt of the manifold 



h{-K,t) hix,t) + ^ 



(5) 



It was observedE3 that the automata of Ref. |l^ present 
strong dependence of the velocity and of the value of fc 
on the parameters m^. In our A^-model, rotational in- 
variance is manifestly broken as an elastic energy 6A*/12 
generates a piece 



(6) 



b V^h{Vhf 



in the equation of motion. From dimensional analysis,EiEj 
among the orientation-dependent differential operators, 
the KPZ term ^ (V/i)^ is more relevant than the term 
eq. (^. The non-linear KPZ term appears in domain- 
growth models |-without disorder, where it is coupled 
to the velocityEj Adding it to eq. (||) constitutes the 
quenched KPZ equation. Its critical exponents have 
not been computed within the functional renormalization 
group, and it is not even clear how this term is generated 
in the absence of a finite velocity, at depinning.Q 
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FIG. 3: Mean square extension W"^ vs system size for the 
(2 + 1) -dimensional elastic manifold for energies i5ci(A) = 
-A'' for 7 — 1.8 (upper, circles), 7 = 2 (middle, squares) and 
7 = 4 anharmonic energy (lower, triangles); M = L. The 
inset shows the mean local elongation (|A|) for 7 = 2 and 
7 = 4. For 7 = 1.8, {|A|) clearly diverges (not shown). Error 
bars are smaller than the symbol sizes. 
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FIG. 4: Mean square extension vs system size (M = 
L) for the (3 + 1) -dimensional elastic manifold for energies 
Bel (A) = ^^A"*" for 7 = 2 (harmonic energy, upper curve) 
and 7 = 4 (lower curve). The inset shows ^^2 and (^^4 vs 
d. The lines are the one-loop and two-loop results for the 
A^-model. 



> 1. This implies that 



Kardar et a/.lill23 have suggested the force anisotropy 
due to the disorder as a possible generating mecha- 
nism, but none of the automata models present such an 
anisotropy. However, in their construction rules, a met- 
ric constraint is hidden. This implements a strong elastic 
potential, and naturally generates terms like eq. (^. In 
our opinion, the main unresolved theoretical issue is to 
understand how this term reduces to the KPZ term under 
coarse-graining. 

As discussed in Ref. ^, there is an easy way to un- 
derstand why the (1 -I- l)-dimensional A^-model is un- 
stable |to higher-order corrections: A roughness C > 1 
impliea^j that the mean local elongation (| A|) grows with 
the system size at least as L''~^, and thus diverges in the 
thermodynamic limit L, M — > 00. Higher-order terms 
in the elastic energy are thus trivially relevant in one 
dimension. □ 

In 2 + 1 dimensions, we have studied samples up to 
size X M = 512^ x 512. For the A^-model (shown 
in Fig. ^), we find excellent scaling for a roughness 
= 0.753 ± 0.002, which again falls between the re- 
sults of the one-loop and two-loop renormalization group 
calculations. Our result improves by an order of magni- 
tude the precision of the previous estimateJl3 Note that 
our algorithm allows us to know, without invoking uni- 
versality arguments, that the exponent Ca^ ~ 0.753 is 
that of the two-dimensional quenched EW model. 

In Fig. H, we also show our results for the A*-model. 
We find a different value for C, namely Ca«i = 0.45 ± 
0.01. This last value is significantly smaller than ( ^ 0.48 
found in Ref. ^ using a cellular automaton. 

The two-dimensional A^-model is not only unstable 
with respect to stronger elastic potentials, but also to 
any weaker ones. We have studied the A'''-model, for 7 
slightly below 2: Already for 7 = 1.8 (c/Fig. we find 
that the roughness exponent changes drastically, as we 
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obtain clear indications that (^1 
the surface breaks. 

We also studied elastic energies ~ vT+A^ in the 
(2+l)-dimcnsional problem. In one dimension, this func- 
tional form corresponds to the length of the string, and 
the A^-model is often understood as the first term from 
an expansion of the length in powers of A. The square 
root is softer than its first expansion coefficient; we again 
find values of Cyi+A' in excess of one for large system 
sizes. 

As in 1 -I- 1 dimensions, we also find for the (2 -I- 1)- 
dimensional A^-model indications of a (much slower) di 
vergence of the mean local extension, as shown in Fig. 
For the A^-model, it seems to diverge logarithmical 
whereas for the A^-model it saturates already for very 
small systems. Again it appears that only the A^-model 
has a proper thermodynamic limit for the critical mani- 
fold. 

We also performed extensive computations in 3 -f 1 
and in 4 -|- 1 dimensions, where we have reached sample 
sizes of X M = 64^ x 64 and x M = 32^ x 32, 
respectively. For the harmonic model in 3-1-1 dimensions, 
we obtain (^2 = 0.355 ± 0.01 again in good agreement 
with the numerical results of Ref. |l^, while our result 
Ca2 — 0.25 ± 0.02 for the A^-model is in contradiction 
to Ref. |l^ which suggested an exponent ( ^ 0.38 larger 
than for the harmonic case. We checked very carefully 
that our estimate for ^ is independent of the parameters 
a and b. 

At last, mir algorithm is able to investigate 4+1 
dimensions .Bj We have computed manifolds of sizes L — 
8, 16, 32 (i = M). Even though we have only three 
points, we find for the A^-model an inconsistent fit to 
a functional form W'^{L) ^ L^^ for sizes L where it was 
already good in 3 + 1 dimensions. At the upper crit- 
ical dimension, one would expect logarithmic behavior 
for (L) . The most natural explanation of our findings 
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0.355 ± 0.01 


0.25 ± 0.02 



TABLE I: Roughness exponents as a function of dimension d 
for the (d -I- 1) -dimensional driven manifold problem, both for 
the A'^-model (harmonic elastic energy) and the A^-model. 



in 4-1- 1 dimensions is that the A^-model's upper critical 



dimension is 



4. This will still have to be confirmed 



for larger sample sizes. 

For the A^-model, we also do not find a good fit for (, 
but the deviations from the functional form W'^{L) ^ L?'^ 
are less striking than for the A^-modcl. It is unclear to 
us whether the roughness exponent of the A'^-model can 
be higher than for the A^-model. 

In conclusion, we have directly computed roughness 
exponents for d-dimensional manifolds in a (d -I- 1)- 
dimcnsional embedding space for d = 1, . . . ,4. Some of 
our findings nicely fit into the existing theoretical frame- 
work. For example, the upper critical dimension due — 4, 
conjectured on the basis of dimensional power counting, 
is consistent with our data for the A^-model. 

In the harmonic greement with the two-loop cal- 



culation is fair, and the exactness of the one-loop result 
can be excluded (c/ inset of Fig. ^). Our numerical work 
relies on no additional hypothesis, as we compute well- 
defined critical manifolds in finite systems. Even for the 
harmonic case, though, we would like to understand the 
divergence of the local scale, for which we find evidence 
even in c? = 2: slightly sub- harmonic potentials lead to a 
break-up of the critical surface. 

For the first time we have presented a Hamiltonian 
model which reproduces the results of the quenched KPZ 
class. The connection between our A'^'-model and the 
equations of motion in the continuum limit is transpar- 
ent. The analysis of the A''-model gives a very precise 
estimate for the value of ^ as a function of dimension d. 
It would be very interesting to understand how the lo- 
cal differential operators, as W^h{Wh)^, which correspond 
to elastic Hamiltonians, renormalize into more relevant 
terms, as the KPZ term. 
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